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Abstract 

The semileptonic decay of a 6-quark, b — > civ, is considered in the relativis- 
tic limit where the decay products are approximately collinear. Analytic 
results for the double differential lepton energy distributions are given for 
finite charm-quark mass. Their use for the fast simulation of isolated lepton 
backgrounds from heavy quark decays is discussed. 
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Many new physics signals at hadron colliders involve isolated hard leptons as a distinguishing 
feature. Sequential decays of supersymmetric particles |I| and Higgs boson decays to Z, W or r 
pairs are but two examples. In all these cases the production of heavy quarks, in particular 
bottom and charm, and their subsequent semileptonic decay constitutes an important back- 
ground. Even though lepton isolation, the requirement that little hadronic energy is deposited 
in the vicinity of the charged decay lepton, can reduce these heavy quark backgrounds by large 
factors, the sheer size of the bb or cc production cross section makes heavy flavor backgrounds 
dangerous [pf. 

For the simulation of such heavy flavor backgrounds the large suppression factors due to 
lepton isolation pose a special problem: large Monte Carlo samples must be generated in order 
to analyze the phase space distributions of the surviving events. While the full five- dimensional 
distribution of & — ^ civ decay is easily implemented in a Monte Carlo program ||, this procedure 
does not always generate isolated lepton events in a sufficiently fast and efficient manner. 

In this brief note I describe how a fast short-cut is provided by analytic expressions for 
the lepton energy distributions in the laboratory frame. Leptons of sufficiently high transverse 
momentum can only result from the decay of very energetic b or c quarks. In turn, this implies 
that the parent quarks must be moving relativistically in the lab, which results in the decay 
products moving approximately collinear to the parent quark direction. In this relativistic limit, 
only the energy fractions of the decay particles, as compared to the heavy quark energy, are 
needed for a full description. 

To be definite, consider the decay b — -> civ and denote the energy fractions of the neutrino, 

the charged lepton and the c-quark by 

E v Ei E c 

&b &b &b 

respectively. Obviously they obey the constraint x + y + z = 1. The smallest energy for the 

charm quark is reached when, in the b rest frame, it is emitted opposite to the 6-quark direction, 

recoiling against a collinear lepton-neutrino pair: 



In the spectator quark model, and for unpolarized 6-quarks, the double differential 6-decay dis- 
tribution can be determined analytically. I find 

1 d 2 T 2c / . , r s . (2 - c) x + c \ 

f^Ty = W)y C{1 - X)lC+i3 - C)X]+3rV 1-x-V J' <3) 

Here 

c= 1 - r - X - y = l-l, (4) 
1 — x — y z 

and /(r) is the phase space suppression factor for the b — > c^z/ decay due to the finite charm 
quark mass |5|f|, 

/( r ) = (1 -r 2 )(l -8r + r 2 ) - 12r 2 logr , (5) 

which is quite sizable for 6-decay: f(r) = 0.42 for r = 0.12. 

For c — > decay the (V — A) x (V — A) structure of the weak decay amplitude implies 
a double differential decay distribution identical to Eq. @, but with the role of charged lepton 
and neutrino energy fractions interchanged, i.e. for charm decay x = Ei/E c , y = E v /E c , z = 
E s /E c >r = m 2 s /m 2 c . 

In the massless limit, m c = i.e. r = and c = 1, the double differential distribution of 
Eq. (H) reduces to 

CT = 2(1 - I)(1 + 2X) ' (6) 

which leads e.g. to the well known lepton decay distribution in r — > iv(V T decay ||[7j 
1 dT /■!-» , 1 d 2 T !,„ /r , , 



r 0<i , a d ^5^=3( 1 -^H 5 + ^-V). (7) 

The double differential decay distribution of Eq. (|3|) provides an adequate description of 
charged lepton and missing transverse momentum distributions and their correlations in typical 
collider physics applications. Because of its simple algebraic form, it may be folded analytically 
with algebraic fragmentation functions, like the Peterson fragmentation function M for 6-quarks. 
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Also, the collinear limit is sufficiently simple to cast phase space limits into limits on the mo- 
mentum fractions of the 6-decay products, once the momentum of the parent b is known, e.g. in 
a Monte Carlo program. 

Another application is the effect of lepton isolation on the observable charged lepton or 
missing transverse momentum distributions in b decays. A typical lepton isolation cut limits 
the energy fraction carried by the charm quark, to e.g. 10 % of the observable lepton energy, 
or imposes an upper limit, e.g. 5 GeV, on its transverse energy. Taking m& = 5.28 GeV and 
m c = 1.87 GeV, i.e. using the lightest meson masses in order to approximately obtain the correct 
kinematics for the heavy quark decays, one finds z > r = 12.5%, which, at face value, excludes 
any events where the charm quark would carry as little as 10% of the charged lepton energy 
or would limit the 6-quark E?, and thereby the maximum lepton E? to about 40 GeV when 
Etc < 5 GeV is required. However, these limits are imposed in the experiment on observed 
hadrons, or calorimeter response in some cone around the lepton direction. For the soft hadronic 
depositions inside the lepton isolation cone, non-perturbative corrections (from fragmentation or 
underlying event contributions) or fluctuations in the calorimeter response lead to considerable 
uncertainties in the true energy fraction z carried by the charm quark. The low energy tails of 
the calorimeter response to charm quarks are largely responsible for fake isolated lepton events. 
As a result, the actual ^-distribution of the charm quark requires detailed simulations, except 
for the general statement that small values of z, close to their kinematic limit z — r, are strongly 
favored. 

The double differential decay distribution derived above allows to assess the effects that z- 
smearing has on the observed lepton distributions. At fixed z = 1 — x — y we may study the 
charged lepton energy distribution 



or the analogous neutrino energy distribution 1/T Z dT z /dx. Here N(z) is a normalization factor 
which is obtained by direct integration of Eq. (^|): 
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FIG. 1. Normalized energy distributions of (a) the charged lepton and (b) the neutrino in 
b —> vie decays in the spectator model. The individual curves correspond to fixed charm quark 
energy fractions z = 0.14 (red), 0.16 (blue), 0.18 (green), 0.2 (magenta), which are slightly above 
threshold, given by z = r = m 2 c /m1 = 0.12. 

= (1 - - )(1 -z) ((1 - rf + \ (1 - -)(1 - z) (4* - r - 1 + 4-)) . (9) 
j(r) z V o z z J 

These neutrino and charged lepton energy distributions are shown in Fig. [I]. They change very 
little with z, in the z-range leading to isolated leptons. The largest z-dependence is found near the 
kinematic limits, somewhat affecting the hardest charged leptons and the softest neutrinos. The 
modest ^-dependence of the lepton distributions implies that the precise ^-distribution produced 
by the lepton isolation cuts is not needed for an adequate description of lepton momentum 
distributions. 

One thus finds that in 6-quark decays which lead to isolated leptons, the charged lepton and 
missing transverse momentum distributions, and their correlations, can be modeled quite reliably, 
making use of the double differential decay distribution described here. A first application in 
collider phenomenology appears in Ref. where bb+ jets backgrounds to H — > rr searches are 
discussed: the decay distributions in the collinear approximation allow to considerably improve 
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Monte Carlo statistics. Similar improvements are foreseen in the simulation of heavy quark 
backgrounds for many new physics signals involving isolated charged leptons. 
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